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ABSTRACT 

We model the population of AM CVn systems in the Galaxy and discuss the detectability 
of these systems with optical, X-ray and gravitational wave detectors. We concentrate on 
the short period (P < 1500 s) systems, some of which are expected to be in a phase of 
direct impact accretion. Using a self-consistent model for the star formation history and radial 
distribution of stars in the Galaxy plus simple models for the emission of optical and X- 
ray radiation from the AM CVn systems and interstellar absorption, we derive the sample 
of short-period AM CVn systems that can be detected in the optical and/or X-ray bands. At 
the shortest periods the detectable systems are all X-ray sources, some with periods as short 
as three minutes. At periods above 10 minutes most detectable systems are optical sources. 
About one third of the X-ray sources are also detectable in the optical band. We also calculate 
the gravitational wave signal of the short-period AM CVn systems. We find that potentially 
several thousand AM CVn systems can be resolved by the gravitational wave detector LISA, 
comparable to the expected number of detached double white dwarfs that can be resolved. We 
estimate that several hundreds of the AM CVn systems resolvable by LISA, are also detectable 
in the optical and/or X-ray bands. 

Key words: gravitational waves - binaries: close - white dwarfs 



1 INTRODUCTION 

AM CVn systems are binaries, with periods less than 1 hour, 
in which a white dwarf accretes matter from a low-mass, 
helium-rich, object. T he mass transfer is driven by gravitational 
wave radiation (e.g. Paczynskj] Il967t) . They are observed at 
optical wavelengths as faint blue and variable objects . For a 
review s of observ a tional work on AM CVn systems, see IWarnetl 
119951) : Solheim 1 2003). A discussion of the formation and 
modelling of the population c an be f ound in lTutuko v &_Yun gelsonl 
ll996t). iNelemans et all fcOOlat). h e reafter INPVYI and 
Podsiadlowski, Han & Rapp aportl 120031). Iffils & Bended 12000 ) 
and Nelemans, Yungelson & Portegies Zwart (2001b) discussed 
the contribution of the AM CVn population to the Galactic 
unresolved gravitational wave background. 

The number of known or suspected systems is rapidly increas- 
ing: six new (candidate) AM CVn systems have been f ound in the 
last fi ve years: CE 315, with a period of 65.1 rmnutejsjRuiz^ydJ 
| 200ll): KL Dra with a period of 25 minutes Jjha et all 1 199a 
IWood et all2002l) : ' '2003aw", with a possible period of 34 minutes 
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IWoudt & Warned" <2003bl): RX J1914.4+2456 (V407 Vul) with a 
period of 9.5 minut es Icropper et alj|l99 8l: iRamsav et aljEoOOt 
IRamsav et all2002tt iMarsh & Steeghsll2002l). KUV 01584-0 939 
(ES Cet) with a period of 10.3 minutes iWarner & WoudJ2002h and 
RX J0806.3+1527 w ith a period of 5.3 minutes ilsrael et all20ol 
IRamsav et al1l2002j) . RX J1914.4+2456 and RX J0806.3+1527 
were discovered as ROSAT X-ray sources iMotch et alJ Il996t 
Ilsrael et all 19991) . For these two sys tems the interpretation as AM 
CVn systems has been questioned <Wu et aljEool iNorton etafl 
2002), but not for the other four new systems. 



The last three of the new systems mentioned above are 
particularly interesting as theoretical models of the evolution of 
AM CVn syste ms and their curre nt Galactic population (e.g. 
iTutukov & Yungelsonlll996llNPVYb show that the expected num- 
ber of systems at periods shorter that ~10 minutes is small, and 
an optically selected sample will be dominated by longer period 
systems. However, at these short periods the mass transfer rate is 
as high as 10~ 6 — 10 -8 Mq yr , and therefore systems might 
be discovered as X-ray, rather than optical sources, as is the case 
for the two ROSAT sources. At these short periods the systems are 
also promising sources for low-frequency gravitational wave (GW) 
detectors (e.g. Israel et al. 2002 ; Marsh & Steeghs 2002). 
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In this paper we determine the range of orbital periods in 
which AM CVn systems could be discovered as X-ray sources, 
using simple, average temperature and blackbody emission mod- 
els and taking into account interstellar absorption. We estimate the 
number of such systems in the Galaxy and their brightness in the 
optical band. We compute the number of short-period (P < 1500 s) 
AM CVn systems that could be detected by the planned low- 
frequency gravitational wave radiation (GWR) detector in space, 
LISA, and focus on the subset of these systems that are expected 
to have optical and/or X-ray counterparts, as combined observa- 
tions will lead to a better understanding of the physical state, for- 
mation and population of AM CVn systems. We will use the term 
'AM CVn system" throughout the paper for all double-degenerate 
interacting binaries with helium secondaries, irrespective of their 
optical properties. 

In Sect. [2] we review our Galactic model of A M CVn sys- 
tems, discuss the changes we made with respect to iNPVYl and 
the uncertainties in this model. The emission of optical, X-ray and 
gravitational wave signals from AM CVn systems is discussed in 
Sect. [3] The selection effects and the resulting expected popula- 
tion of sources detectable in the optical and X-ray bands, and as 
GWR sources are presented in Sect. [4] In Sect.|5]we consider the 
LISA sources that have optical and/or X-ray counterparts. Discus- 
sion (Sect.|SJ and conclusions (Sect.Q follow. 
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Figure 1. The distributio n of stars in the x — z pla ne of the Galaxy for our 
implementation of the Boissier & Prantzos 1 1999) model at Galactic ages 
of 1, 6 and 13.5 Gyr, using a number of tracer particles proportional to the 
total mass in stars (M„ ) in the Galaxy as indicated in the plot. The change 
in the radial extent of the star forming region in time is clearly visible. For 
the total SFR as function of time, see Fig.l2l 



2 A MODEL FOR THE CURRENT GALACTIC 
POPULATION OF AM CVn SYSTEMS 

2.1 Formation of AM CVn systems 

For a detailed discussion of the formation of AM CVn systems we 
refer to NPVY, and references therein, so here we only give a brief 
summary. There are two main routes for the formation of AM CVn 
systems: (i) via a phase in which a double white dwarf loses angu- 
lar momentum due to gravitational wave radiation and evolves to 
shorter and shorter periods to start mass transfer at periods of a few 
minutes and (ii) via a phase in which a non-degenerate helium star 
transfers matt er to a white dwarf ac c retor. 

Recently Podsiad lowski et all 120031) proposed a third sce- 
nario to form AM CVn binaries, similar to the usual scenario for 
the formation of cataclysmic variables (CVs), via a phase of mass 
transfer in the course of which a low-mass star with initially al- 
most exhausted hydrogen core or a small helium core becomes a 
helium white dwarf 1 . We found that with the assumptions used in 
this paper (see below) this scenario is relatively unimportant for the 
short-period systems discussed here (giving less than two per cent 
of the systems), so for the remainder of the paper we leave it out of 
our model, but briefly come back to it in the discussion (Sect.|6j. 

One of the main conclusions of NPVY was that the uncertain- 
ties in the models for the formation of AM CVn systems resulted 
in widely varying numbers of expected syste ms in the Galaxy. The 
most optimistic and most pessimistic cases of iNP VYl gave roughly 
80 and 18 million systems in the Galaxy at present. The uncer- 
tainty is particularly severe for the shortest period AM CVn sys- 
tems, which probably have to form through the double white dwarf 
scenario. In lNPVYl we showed that at the start of the mass trans- 
fer between double white dwarfs there is a phase of "direct impact 
accretion", in which the gas stream flowing from the donor to the 

Such a scenario has been proposed before (e.g. iTutukov et all 19871) as a 
formation scenario for ultra-compact X-ray binaries 



accretor hits the accr etors surface di rectly without forming an ac- 
cretion disc (see also lWebbinkl 19841) . In this phase the stability of 
the mass transfer depends crucially on the tidal (or other) c oupling 
between the spin of the accretor and its companion JnpvyI see for 
a deta iled discussion of the stability [Marsh. Nelemans & Steeghsl 
2003). 

The formation of systems formed via the helium star scenario 
is also uncertain because of the possibility of so called edge-lit- 
detonation (ELD) of the accreti ng white dwarf before the system 
reaches the AM CVn phase (see lNPVYl) . I n the cur rent study we 
use the most optimistic model (Model II of INPVYl) . in which we 
assume efficient coupling in the direct impact phase, and that accre- 
tion of at least 0.3 Mq of He is necessary for ELD. But we note that 
recent calculations, including rotation, of helium accretin g white 
dwarf s suggest that ELD might not happen at all lYoon & LangeJ 
2003). This would increase the formation rate of AM CVn systems 
via the helium star scenario in our most optimistic model by about 
15 per cent. For the short period systems discussed here the total 
uncertainty is probably of the order of a factor 10 at least. 

We note that we do not take into account the unstable helium 
burning on the surface of carbon-oxygen white dwarfs accreting 
helium at high rates and the possible associated mass loss (see, e.g. 
lKato&Hachisull999h . 



2.2 Galactic model, star formation history and stellar 
evolution 

Instead of the simple Galactic model we used before, we imple- 
mented a star formation rate (SFR) as function of time and position 
in the Galaxy based on the model of Bois sier & PrantzosM 19991) . 
We obtain the SFR as function of R and t by linear interpolation 
between values in a table (S. Boissier, private communication), and 
assume that the probability of a star being born with a radius R or 
smaller goes with the integrated SFR, i.e. 
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Figure 2. Illustration of the difference in Galactic model between the current work (left panel) andlNPVYi (right panel). Shown are the SFR (solid lines), 
the birth rate of AM CVn systems (dashed line, left ordinate), and their total number (dot-dashed line, right ordinate) as function of time in the Galaxy. 
The differences in the birth rate and number of AM CVn systems reflect the differences in star formation history and the treatment of stellar evolution (see 
Sect. 12.3^ Note that in the exponential model the Galactic age is 10 Gyr, while the Boissier & Prantzos 1 1999) model has an age of 13.5 Gyr. 



P(R,t) 



f*SFR(R',t)2TrR' dR' 
j-flma* SFR (#, t ) 2 7T R> dR' ' 



(1) 



where 7? max is the maximum extent of the disc (19 kpc) and we 
assume that the stars do not migrate radially in time. The age of 
the Galaxy in this model is 13.5 Gyr. As the Boissier ~& Prantzosl 



( 1999) model gives the SFR projected onto the plane of the Galaxy, 
we assume a ^-distribution, 

P{z) <x sech(z/2 h ) 2 , (2) 

as in lNelemans et alJfcOOld) . where zu = 200 pc, neglecting its age 
and mass dependence. We have added a bulge to the model, simply 
by doubling the SFR in the inner 3 kpc of the Galaxy compared to 
iBoissier & Prantzosl fT999). As a result the total mass in this region 
is 2.6 10 10 Mq at the current age of the Galaxy, co nsistent with 
kinematic and micro-lensing results (e.g. iKlypin et alj|2002l) . We 
distribute the stars in the bulge according to 



Pbuigc = exp(-(r/0.5kpc) ), 



(3) 



where r = \J x 2 + y 2 + z 2 . 

In Fig.Q] we plot the resulting distribution of stars in the x — 
z plane of the Galaxy, showing the ^-distribution of stars and the 
change in the radial extent of the star formation as function of time. 
For the SFR integrated over R as function of time, see Fig.0(left 
panel). 

The important effect resulting from the different choice of 
Galactic model is that f or most of the time in the G alactic history, 
the SFR is lower in the B oissier & Prantzosl 1 19991) model than in 
the 'exponential' model used in lNPVYl (see Fig.|2|; in particular, 
the current SFR is lower than the one we obtained previously (1.9 
vs 3.6 Mq yr" 1 ). This has a strong effect on the current birth rates 
of AM CVn systems, since the latter follow star formation history 
with a typical delay o f several Gyr iTutukov & Yungelsonlll996t 
iNelemans et al 12 001 a). However, the integrated SFR remains quite 
similar (8.2 x 10 10 vs. 8 x 10 10 M ), partly due to the different 
Galactic ages (13.5 vs. 10 Gyr). Another characteristic of our new 
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Figure 3. Histogram of distances to the Earth of short period AM CVn 
systems for our implementation of the IBoissier & Prantzosl fl 999 ) model 
with added bulge and for an exponential radial density distribution. 



model is that the Galactic distribution is more centrally concen- 
trated than with an exponential radial density distribution , mainly 
due to the inclusion of the bulge. In Fig. [5] we show the distribu- 
tion of distances from the Earth of the population of short-period 
AM CVn systems in the Galaxy for the two models. The concen- 
tration in the Galactic centre with the new model is clearly visible. 
As the longer period AM CVn systems are older, their distribution 
is even more centrally concentrated 

Another new element in our current study is th e use of the 
IM F as pro posed by Krouoa, Tout & Gilmore ( 1993), rather than 
thelMiller&ScalcMl979l) IMF as used in lNPVYl Because the 
lKroup"a^^jT^993l) IMF is steeper for masses above 1 Mq , the 
same star formation history leads to fewer AM CVn systems, es- 
pecially for the helium star formation scenario, the efficiency of 
which is reduced by about 25 per cent (see Sect. 12. 3> . 
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(lO" 3 yr" 1 ) 


(10 7 ) 


double white dwarfs 


1.3 


2.3 


helium stars 


0.27 


1.1 


total 


1.6 


3.4 



Table 1. Present day birth rates and total number of AM CVn systems in 
the Galaxy in our current model, for the different formation channels. 



We use the population synt hesis code SeBa, d eveloped b y 
Porte eies Zwart & VerbuntHl996l see also lNelemans et all2001ch . 
We assume a 50 per cent binary fraction, an initial distribution of 
separations in binary systems (a) flat in log a between contact and 
10 6 Rq, a flat mass ratio distribution and a thermal eccentricity 
distribution. 2 The current study is based on a model with 250,000 
ZAMS binaries. For the treatment of the stellar evolution and the 
mass transfer we refer to our earlier papers, but recall that we use a 
common envelope formalism based on the angular momentum bal- 
ance for the first episode of unstable mass transfer in the system, as 
detailed in lNelemans et alj )2000. 2001c). Two small changes that 
turn out to be important are in the formation and evolution of white 
dwarfs. We updated our treatment of the growth of the core in gi- 
ants as described in Nelemans et al. (2001c), by including a simple 
model of the second dredge-up to be more in line with detailed stel- 
lar evolution models (e.g. Hurlev et al. 2000). The second change 
is in the treatment of white dwarf cooling (see Appendix 1X1 and 
white dwarf mass - radius relatio n, for which we n ow use the fit by 
Eggle ton, as given in |Vgrbjint_&.Rappaport (1988), rather than the 
Zanols kv&Salpeterl ri969) relation, which is not appropriate for 
massive white dwarfs. 

The last change we made compared to lNPVYl is in the max- 
imum rate at which a main sequence star can accrete before it ex- 
pands and mass will be lost from the system. We used to limit this 
to accretion at the thermal time scale o f the accre t or. Ho wever, as is 
shown e.g. by Kippen hahn & Mever-Hofmeistedll977l) . higher ac- 
cretion rates do not necessary lead to expansion of the accretor. We 
therefore use the treatment of Pols & Marinus 1 1994) which takes 



the reaction of the radius of the accreting star upon mass gain into 
account. 



2.3 Total population of AM CVn systems in the Galaxy 

With the changes described in the previous section, we obtain, for 
the two formation channels discussed here, the current birth rates 
and total number of AM CVn system s in the Galaxy as shown in 
TableQand Fig.[2](compare Table 1 in lNPVY}). Th e difference be- 
tween the numbers presented here and in InpvyI are due to the 
use of the new Galactic mode, different IMF and the changes in 
the treatment of stellar evolution described in Sect. 12.21 All these 
change the birth rate of AM CVn systems at any time, e.g. at the 
current Galactic age the birth rate of the white dwarf family goes 
down by ~60 per cent due to the Galactic model and is reduced 
further by another ~40 per cent due to the changes in treatment of 
white dwarf formation and mass - radius relation. The birth rate of 



P(s) 



Figure 4. Histogram of periods of the shortest period AM CVn systems. 
The shaded histogram gives the number of systems which are in the direct 
impact phase (but note the logarithmic y-scale). The observed systems are 
indicated at their suggested orbital periods. The number of systems in these 
bins are indicated. 



the helium star family goes down by ~60 per cent due to the Galac- 
tic model, ~25 per cent due to the IMF and ~40 per cent due to 
the higher accretion rate limit for main-sequence stars (because less 
systems pass through common envelopes and become close enough 
to produce AM CVn systems). Partly compensated by the change 
in the age of the Galaxy this results in 55 and 65 per cent fewer 
systems present in the Galaxy for the two families, respectively. 

In Fig.|4]we show a histogram of the shortest period systems 
as function of orbital period, indicating the systems in the direct 
impact phase and the periods of the observed (candidate) AM CVn 
systems. 



3 MODELLING THE OPTICAL, X-RAY AND GW 
EMISSION OF AM CVn SYSTEMS 

In order to be able to compare our results with (future) observa- 
tions we have to simulate the emission produced by our model 
systems. We determine the temperature and luminosity of the dif- 
ferent emission sources and determine how much energy is emit- 
ted in a certain band, assuming the emitted spectrum is a black- 
body. For the present study we only discuss emission in V-band 
and the ROSAT 0.1-2.4 keV X-ray band plus gravitational wave 
radiation. We chose the ROSAT band, because the two candidate 
systems were found as ROSAT sources and because we can com- 
pare our results to the ROSAT all-sky survey. However, in Sect. [6] 
we briefly discuss the expected results for the Chandra and XMM. 

There are four main emission sites: the accretion disc and 
boundary layer between disc and accreting white dwarf, the im- 
pact spot in the case of direct impact accretion, the accreting star 
and the donor star. 



2 All binaries discussed here will be circularized to somewhat smaller sep- 
aration before the first phase of mass transfer and thus the shape of the initial 
eccentricity distribution hardly influences the number of AM CVn systems. 



3.1 Optical emission 

The disc. We assume that half of the accretion luminosity is ra- 
diated by the disc (the other half in the boundary layer), i.e. 
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Figure 5. Absolute magnitudes and temperatures of the different optical 
emission components as function of orbital period. The limits of the regions 
are defined by systems with initial masses of 0.15 + 0.25 Mq and 0.3 + 1.0 
Mq . For the accretor and donor we give the upper limits on the magnitude 
and temperature for the extreme case that the donor/accretor white dwarf 
was young and hot at the start of the AM CVn phase, i.e. at the top of the 
cooling curve (see Appendix lAl. 



Ldis 



1 

~R^i 



erg s 



(4) 



where M and R are the mass and radius of the accretor and Rli 
is the distance of the first Lagrangian point to the centre of the 
accretor and m is the mass transfer rate. Because the two white 
dwarfs are so close together the usual assumption that the matter 
falls in from infinity is not valid ( Han & Webbink 1999). We model 
the op tical emission from the disc like in lNPVYl according to lWadd 
( 1984), as a single temperature disc that extends to 70 per cent of 
the Roche lobe of the accretor, radiates as a blackbody and use the 
bolometric correction appropriate for its temperature. 

The donor. We model the emission from the donor as the cooling 
lu minosity of the white dwarf, using approximations to the models 
of H anse nl ll99i . see the Appendix IAI 

The accretor. We model the emission from the accretor as being 
the unperturbed cooling luminosity of the white dwarf. 

In our simple model we do not take into account such effects 
as irradiation of the donor by the impact spot or disc, he ating of 
the accretor by radiation from the disc/boundary layer (e.g. Pringle 



1988), compressional heating (e.g. lTownslev & Bildsten l 2002l) or 
tidal heating of the two white dwarfs prior to contact llben et alj 
Il998h . 

In Fig.|5|we show the range of absolute V-magnitudes (Mv) 
and temperatures of the different optical components described 
above, for systems forming from double white dwarfs. The top 
panel shows the range of absolute magnitudes of the discs as func- 
tion of period (in grey). The two symbols show the observed abso- 
lute magnitudes for AM CVn and GP Com [using distances of 235 
pc for A M CVn (C. Dahn 3 , p rivate communication) and 70 pc for 
GP Com (Thorstensen 2003)], while the bars in the bottom panel 
show the range of effective temperatures for these systems. The up- 
per and lower limit of the shaded regions are for two combinations 
of masses that bracket the range of masses of the majority of the 
systems: 0.3 +1.0 and 0.15 + 0.25 Mq respectively. 

Because for the accretor and donor we assume that their lumi- 
nosity is the cooling luminosity, their magnitudes depend on their 
ages, rather than on the orbital period. However, during the AM 
CVn phase there is a direct relation between time and period (e.g. 
figure 4 in lNPVYl) so for given cooling luminosity at the start of 
the mass transfer we can plot the magnitude of the cooling white 
dwarf as function of orbital period. In Fig.|5|we plot that curve for 
the extreme case in which the cooling luminosity at the starts of 
the mass transfer is that of a newly born, very hot and luminous 
white dwarf. In reality the donor and accretor are typically of the 
order of a Gyr old, so they evolve along a track that lies below these 
curves. Towards the longer periods, the age of the system is a few 
Gyr anyway, so the curves for all initial luminosities converge. 

The few observed systems fall in the grey area in Fig.|5|indi- 
cating that our simple model is at least able to understand these sys- 
tems satisfactorily except the temperature of GP Com 4 The same 
model can be adopted for the systems that form via the helium-star 
donors. In these cases, however, the curves shift to longer periods 
(see lNPVYl) . 



3.2 X-ray emission 

In modelling the X-ray emission of AM CVn systems we distin- 
guish two cases: with and without disc. 

If there is no disc, the gas stream hits the companion directly 
and a small area of the surface will be heated. We assume that the 
full accretion luminosity is radiated as a blackbody with a temper- 
ature given by 



Tbb 



7- 



(5) 



where L acc and R are in solar units and L acc is defined by 
Eq. 0. The fraction / of the surface that is radiating depends on 
the details of the accretion, but will in general be small. We used 
/ = 0.001, which is consistent with the observed X-ray emission 
of V407 Vul and simple arguments regarding the impact stream 
iMarsh & Steeehi2002l) . 

If there is an accretion disc, the X-ray emission can come 
from the boundary layer between the inner disc and the accretor, 
where we assume half of the ac cretion luminosity is radiated, as 
suggested e.g. by IPringld ll977l) . or it might have a non-thermal 



3 On behalf of the USNO CCD Parallax Team. 

4 Indeed the continuum emission of GP Com originates at the reheated 
accreting white dwarf, rather than in the disc (see Bildsten et al„ in prepa- 
ration). 
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plot, where the evolution is so fast that hardly any systems are ex- 
pected to be observed. 

At periods larger than 1000 - 2000 s the mass transfer rates 
drop to below 10~ 10 Mq yr _1 , the temperature drops below 10 5 K 
and hardly any X-ray emission is produced in the boundary layer. 
The inferred temperatures and X-ray luminosities of RX J0806 and 
V407 Vul are plotted as the bars in the figures. For RX J0806, we 
use the X-ray flux (4 x 10~ 10 erg cm" 2 s -1 ), distance ( 0. 1 - 1 kpc) 
and blackbody temperature (6^4 3 x 10 5 K) as given bv llsrael et alJ 
(2002). For V 407 Vul blackbody temp eratures of 43 and 56 eV 
are derived by Motch et al. 1 1996) and Ra msav et alJ l2002bl) re- 
spectively, while the latter authors derive an X-ray luminosity of 
1.21q 3 x 10 33 erg s _1 , for a distance o f 100 pc and we plot these 
limits for a distance range 100 - 400 pc l Ramsayet al j200ol> . The 
X-ray luminosity of AM CVn is taken from lUlhUl995h . Except 
for the luminosity of RX J0806, our very simple model agrees well 
with the observations. 



3.3 Gravitational wave emission 

We assume that AM CVn systems emit GWR as two point sources, 
i.e. neglecting any effect of the deformation of the donor, which 
fills its Roche lobe and the mass in the accretion disc, if present 
(see R ezzolla et all20oH for a discussion of the small effect of de- 
formation for cataclysmic variables). The amplitude of the GWR 
signal, as can be detected by LISA, is given by (e.g. lEvans et"all 

EH 

where M = (M m) 3/5 /(M + m) 1/5 is the chirp mass, P rb is the 
orbital period and d is the distance to the system. 



Figure 6. Total accretion luminosity and luminosity and temperatures of 
the X-ray emission components (impact spot or boundary layer) as function 
of orbital period. The limits of the regions are again denned by systems with 
initial masses of 0.15 + 0.25 Mq and 0.3 + 1.0 M Q . 



origin. The average temperature of the disc (Wa dell9'84l) . which we 
use to estimate the optical emission is never high enough to produce 
X-rays <Bathetalll974l) . Even though from observations of cata- 
clysmic variables it is clear that the real X-ray emission mechanism 
is mu c h more complicated (e.g.lHoare & Drewll9 91 : Mauch^UtlJ 
Il99lt Ivan Teeseling & Verbuntlll994t Ivan Teeseling et alJll996h . 
we model the X-ray emission of systems with a disc for t he mo- 
ment as coming from the boundary layer, with temperature I Pringle 
119771) 



T BL = 5x10' 



M 

Me 



\W 18 g/s / 

In Fig. [6] top panel, we show how the luminosity in X-rays 
compares to the total accretion luminosity. For the X-ray luminos- 
ity we consider here the emission in the ROSAT band, from 0. 1 - 
2.4 keV. The temperature of the impact spot or boundary layer is 
shown in the bottom panel. The sharp drop in the X-ray luminosity 
and temperature at the bottom of the darker region at log P ~ 2.9 
is caused by the change from direct impact to disc accretion in the 
low-mass system. We assume that the transition from the impact 
temperature to the boundary layer temperature is instantaneous. 
This transition for the more massive system occurs just left off the 



4 RESULTS I: THE DETECTABLE POPULATION OF 
SHORT-PERIOD AM CVn SYSTEMS 

4.1 Selection effects 

To compare our models with the number of observed systems, we 
have to decide which of the model systems would be detectable. 

For the optical emission we focus on model systems that are 
brighter than V = 20 mag, which is similar to the optical mag- 
nitude of the observed short period AM CVn candida tes. To take 
the interstellar reddening into account we use the Sandage <1972l) 
model 



A v (oo) = 0.165 
= 



(tan(50°) -tan(fc)) 
sin(b) 



for b < 50° 

for b > 50°, (8) 



and Av(d) follows from our Galactic model [Eq. J2j] 



Ay(d) — Av(oo) tanh 



dsin(6) 



2/, 



(9) 



For the sources detectable in X-rays we select the systems 
with an X-ray flux in the 0.1 - 2.4 keV (ROSAT) band higher than 
10 -13 erg s _1 cm -2 . To estimate this flux we use the intrinsic flux 
in this band (Sect. l3~2l . the dist ance and an estimate of the Galac- 
tic hydrogen absorption (Morrison & McCammon 1983). We mod- 
elled the hydrogen column density after IPred ehl & Sc hmittl 1 19951) 
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Figure 7. Period vs distance distribution of systems detectable in the X- 
ray and/or optical bands. Top panel: Systems detectable only in X-rays 
(pluses), and systems that are detectable in the optical and X-ray band 
(filled circles: direct impact systems with a bright donor, filled triangles: 
systems with an optically bright disc, filled squares: systems with bright 
disc + donor). Bottom panel: Systems only detectable in the optical band 
(circles: direct impact systems with a bright donor, triangles: bright discs, 
squares: bright disc + donor). 



Nh = 0.179 A v (d) 10 cm" 2 . (10) 



4.2 Systems detectable in the optical and X-ray bands 

In Fig. we show the sample of the systems that would be de- 
tectable in the optical and/or X-ray bands, taking into account the 
selection effects described above. The relevant numbers are given 
below in the text, between square brackets (see also Table[3}. (We 
should recall that the results we discuss represent one possible 
random realization of the model for current Galactic population 
of AM CVn systems, so all numbers given are subject to Pois- 
son noise). We limited ourselves to the short period systems, with 
P < 1500 s. At longer periods the optical emission from the disc 
will completely dominate and the detectable sample has been dis- 
cussed in NPVY. Fig.0 top panel, shows the systems detectable in 
X-rays. Pluses are [220] systems only detectable in X-rays, while 
the filled symbols mark systems also detectable in the l/-band. 

At the shortest periods a significant fraction of the small num- 
ber of systems that exist (see Fig. [4) produces enough X-rays that 
they are detectable even though the sources are close to the Galac- 
tic centre (Fig.Q. At longer periods the X-rays get weaker (and 
softer) so only the systems close to the Earth can be detected. Sev- 
eral of these systems are also detectable in the optical band (filled 



symbols). The [3] filled circles show the systems in the direct im- 
pact phase, where the optical emission comes solely from the donor 
as the impact spot is too hot to emit in the V-band. Systems with 
optical emission from a disc are shown as triangles [75]. There are 
some [28] systems that are close enough to the Earth that the donor 
stars can be seen as well as the discs (filled squares). Above 600 s 
the systems descending from helium star donors appear and have 
high enough mass transfer rate to be X-ray sources, the closer ones 
also being visible in the V-band (filled triangles/squares, as these 
systems always have a disc). 

The X-ray sources at larger distances (pluses) are almost all in 
the bulge and thus are concentrated in a small area on the sky. The 
closer systems, which often are also detectable in the V-band are 
more evenly distributed over the sky. 

The bottom panel of Fig. Q shows the 1230 'classical' 
AM CVn systems, detectable only by the optical emission, for most 
systems only from their accretion disc (open triangles [1057]). Of 
this population the ones closest to the Earth in some cases [169] 
also have a visible donor (open squares). A few [4] direct impact 
systems only have a visible donor, but are not X-ray sources (open 
circles). The majority of the optically detectable systems repre- 
sented in the lower panel of Fig. (with orbital periods between 
1000 and 1500 s) are expected to show outbursts due to the viscous- 
thermal disc instability fTsugawa & O saki 1997) which could en- 
hance the chance of their discovery. 

Comparing these results with the tentative classification of 
V407 Vul and RX J0806 as direct impact systems, detectable in X- 
rays and (just) in the optical band, shows that a few of such systems 
(filled circles) could be expected in the Galaxy, given all the uncer- 
tainties in population synthesis model, the emission model and the 
selection effects. Furthermore it shows that many more short period 
AM CVn systems might be detectable only as X-ray sources. Some 
have periods as short as three minutes. About one third of the X-ray 
sources is also detectable in the optical band. 



4.3 Gravitational waves from short period AM CVn systems 

In lNelemans et alJfcOOlbl) we presented a study of the gravitational 
wave signal from Galactic binaries in which both stars are compact 
objects. In that study we estimated the contribution of AM CVn 
systems to the double white dwarf noise background and found that 
AM CVn systems hardly contribute, even though at frequencies be- 
tween 0.3 and 1 mHz (periods between 33 and 105 minutes) they 
outnumber the detached double white dwarfs. This is because at 
these frequencies their chirp mass has fallen far below that of a 
typical detached system (see also Hils & Bendet 2000). We there- 
fore did not consider them anymore. However at shorter periods, 
where the number of AM CVn systems is much smaller their chirp 
mass is similar to that of the detached systems from which they 
originate. 

In Fig. [8] we show the Galactic population of detached dou- 
ble white dwarfs and AM CVn systems that can be resolved by 
LISA. There are <^11 000 resolved double white dwarfs as well as 
^11 000 resolved AM CVn systems. Their distribution is shown 
by the grey shading. The 200 strongest sources of each type are 
shown as dots to enhance their visibility. Figure|8|cleaiiy shows that 
AM CVn systems may be important sources for LISA. Compared 
to our previous study the 'confusion limit', i.e. the freq uency at 
which the unresolved background disappears (see Nelemans et al. 
(2001b) for details) moves to slightly higher frequency because 
in our current model there are 213 million detached close double 
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Figure 8. Gravitational wave amplitude h as function of the frequency of the wave for the resolved double white dwarfs (left panel, 10697 objects) and 
AM CVn systems (right panel, 1 1000 objects.). The grey shades give the density distribution of the resolved systems normalised to the maximum density 
in each panel (1548 and 1823 for the double white dwarf and AM CVn panel respectively). The 200 strongest sources in each sample are shown as the dots 
to enhance their visibility. In the AM CVn panel the periods of the observed short period s ystems are indicated by the vertical dotted lines. Solid line shows 
the average background noise due to detached white dwarfs as in Nejemanse^d 12001 tat) . The LISA sensitivities for a integration time of one year and a 
signal-to-noise ratio of 1 and 5 are indicated by the dashed lines lLarson et alfcOOd) . 



white dwarfs, instead of 1 10 million before, due to the use of dif- 
ferent IMF, Galactic model and the changes in stellar evolution. 

The crucial difference between AM CVn systems and de- 
tached double white dwarfs in the short period range is the sign 
of their period derivatives. The periods of AM CVn systems in- 
crease due to mass transfer, while the periods of detached systems 
decrease due to GW losses . Assuming conservative mass transfer, 
the change of the GW frequency (/ = 2/P) of AM CVn systems 
formed from double white dwarfs can be derived from 

T = - p = -77 - = -3— 3(1-5) — , (11) 

j P 2 a J or b m 

where m is the mass of the donor star and q = m/M. We can use 
the equilibrium mass transfer rate as derived in iMarsh et alH2003t) 
to get 

/ = -5.8x 10~7 11/3 ( £-) ( 5^: C2 „ 1 Hz s 

where £2 is the logarithmic derivative of the secondary radius with 
respec t to mass and C tTL — dln(.Ri,/a)/dlnm ~ 1/3 (see 
Marsh et all2003l) . We wrote the equation in the same form as the 
corresponding equation for frequency change in detached systems 
(e.g. Eq. (6) in lNelemans et all2001bl) . with the (positive) term in 
parentheses giving the change with respect to the detached case. 
This term is of order unity at the onset of the mass transfer (i.e. 



5 Note that systems formed from helium stars and evolved secondaries 
in C Vs have decreasing periods before they reach their period minimum 
(e.g. Tutukov & Fedorova 1989; Podsiadlowski et al. 2003). However, in 
this phase they are (much) brighter than after the period minimum and are 
not consistent with being classified as "classical" AM CVn systems. 



at highest frequency, above ~5 mHz), decreasing to about 0. 1 at 2 
mHz. 

LISA will be able to measure the frequency change of a num- 
ber of the resolved binaries. In Nel emans et alj 1200 lbl) we esti- 
mated the number of systems for which this can be done by re- 
quiring that in the course of the experiment the frequency would 
change by one frequency bin. However, it really is the phase shift 
resulting from a freq uency change that will be measured (e.g. 
IWebbink & Han 1998), so the requirement has to be that over the 
course of the experiment the phase shifts by a certain amount, com- 
pared to a system without frequency change. Writing the phase as 
an expansion in /, / and /, this means / can be measured if 

A$ = / ^ > A$ min , (13) 

where A$ m i n is the minimum phase shift that can be measured 
^rid T is the duration of the experiment. Analogously / can be mea- 
sured if 

A$ = /— > A$ mi „. (14) 

Measuring the frequency change provides important information, 
as the chirp mass can be measured directly for the detached sys- 
tems. Combined with the wave amplitude h [Eq. Q], this yields 
the distance. For AM CVn systems this argument does not hold, 
as / is determined by the mass transfer [Eq. J12H . However, for 
an assumed mass-radius relation of the donor star, both the model 
mass transfer rate and the mass of the donor star (and thus the chirp 
mass of the syste m to within a factor of ^ 2) as function of pe- 
riod are fixed (e.g. lVilall97lllFaulkneill97ll) . Thus a measurement 
of / can test mass-radius models, which is especially interesting 
because more realistic mass- radius re lations for low-mass white 
dwarfs become available (e . g . 1b ildstenl2002l) . In Table|2|we show 
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Table 2. Number of systems for which / or / can be measured as function 
of the minimum phase shift needed for detection. The duration of the LISA 
mission is assumed to be 1 year, except for the last row, where a mission 
time of 5 yr is assumed. 

the number of systems for which we can measure the frequency 
change as function of the minimum phase shift. 

Finally it would be very interesting to measure the second 
derivative of the frequency, as that gives a direct check whether 
the change in frequency is due to GWR only, or whether other ef- 
fects (like tidal torques etc.) are important iWe bbink & Hani 19981 
Phinney, private communication). For pure GW losses we have 6 

11 f 2 
3/ " 

As can be seen from Table |2| the number of systems for which 
the second derivative can be measured is small, even for a more 
optimistic assumption of a total LISA operation time of 5 yr. For 
none of the AM CVn systems we expect its second derivatives to 
be measurable. This is due to the fact that the AM CVn systems 
form from pairs of low -mass (M < 0.5M Q ) white dwarfs (e.g. 
Nele mans et alj2001al) which have larger radii and thus do not pen- 
etrate into the highest frequency range where / is large (see Fig.[3}. 
In this range there are only massive (M > 1M Q ) detached double 
white dwarf pairs that merge rather than become AM CVn systems. 



5 RESULTS II: SHORT-PERIOD AM CVn SYSTEMS 
DETECTABLE BY LISA WITH OPTICAL AND/OR 
X-RAY COUNTERPARTS 

In Figure |5| we show the number distributions vs. orbital periods 
for LISA sources that have optical and/or X-ray counterparts, di- 
rect impact systems that are GWR and X-ray sources, and the total 
number of systems that are optical and/or X-ray sources (see also 
Table |3j. For comparison, we indicate also the distribution of the 
total number of AM CVn systems with P < 1500 s as the white 
histogram in each panel. The systems with different properties are 
marked by different hatching styles. We will discuss these groups 
now in more detail. 

Resolved LISA sources (GWR) 
The top left panel shows the 11,000 systems that are expected to 
be resolved by LISA (grey histogram, also shown in Fig. [5}, most 
of which do not have optical or X-ray counterparts. There are 151 
short-period systems that have X-ray counterparts (GWR+X) and 
143 systems, mainly with periods between 500 and 1000 s, that 

6 Note that lWebbink & Hani Il998l> give 5/3 as numerical factor, which is 
appropriate for the corresponding equation for the second derivative of the 
orbital period. 



have optical counterparts (GWR + Opt). Only 3 1 systems have both 
optical and X-ray counterparts (GWR + Opt + X). 

Direct impact accretors 
The top right panel shows the 2680 systems that are in the direct im- 
pact phase of accretion (grey histogram, as in Fig.[4j- Most of these 
systems (2057) are expected to be resolved by LISA [DI (GWR)], 
while at the shortest periods there are 105 X-ray systems that are in 
the direct impact phase [DI (X)]. 

Optical sources 

The bottom left panel shows (in grey) that for periods below 25 
minutes there are 1336 systems that have V < 20 mag (grey his- 
togram). These are the systems shown in the bottom panel of Fig.0 
and are the local population of AM CVn systems. Some (106) of 
these systems are also X-ray systems (Opt + X) and are shown as 
the filled symbols in the top panel of Fig. Q Many of the short 
period systems will show eclips es (even more than the detached 
systems which are discussed in lCoorav et aljEo03h which would 
compliment the GWR data and provide information on the radii of 
the components and the inclinations of the systems. E.g. a crude es- 
timate for a typical system with initial parameters 0.25 + O.6M0 
gives a probability for eclipsing the accreting star of about 30 per 
cent at P — 1000s, and even higher for eclipsing (part of) the ac- 
cretion disc. 

X-ray sources 

Finally, the bottom right panel shows the 326 X-ray sources that 
are also shown in the top panel of Fig.0 The majority is only de- 
tectable in X-rays and many are in the direct impact phase. Only a 
few direct impact systems are detectable in the optical band, by the 
radiation from the cooling white dwarf (the filled circles in Fig.0. 
At longer periods we expect detectable X-ray emission from the 
boundary layers of systems formed from helium stars, some of 
which are also detectable in the V-band (Opt + X). These X-ray 
sources in the ROSAT band in our model are concentrated in the 
Galactic plane, especially towards the Galactic centre. There are 
193 systems (60 per cent of the total) in an area between Galac- 
tic longitudes —50 and +50 degrees and Galactic latitude between 
-10 and +10 degrees. The ROSAT all-sky survey 7 found ~4500 
sources in this region, so the short-period AM CVn systems can 
only account for a small fraction of these sources. 

It is clear from Fig.|5|that LISA is expected to be very effective 
in detecting the short-period AM CVn systems. A large fraction of 
these will be in the direct impact phase. At short periods, a signif- 
icant fraction of the LISA sources might be detectable in X-rays, 
while at longer periods a small fraction might be optical sources. 

LISA sources with optical and/or X-ray counterparts are im- 
portant, as LISA will measure a combination of all the parameters 
that determine the GWR signal (frequency, chirp mass, distance, 
position in the sky and inclination angle, e.g. Hellinas 2003), so if 
some of these parameters (period, position) can be obtained from 
optical or X-ray observations the other parameters can be deter- 
mined with higher accuracy. This is particularly interesting for the 
distances, inclinations and masses of the systems, which are very 
difficult to measure with other methods. 



6 DISCUSSION 

We stress here that AM CVn systems form from binaries with a 
rather small range of initial parameters and thus that their birth 

7 http ://w w w. xray . mpe.mpg . de/rosat/survey/rass-f sc/ 
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Figure 9. Histograms of the population of short-period AM CVn systems, subdivided in different types. In each panel we show the total AM CVn population 
as the white histogram. In the top left panel we show the systems that can be resolved by LISA in grey and subdivide these in the ones that have optical 
counterparts (GWR + Opt), X-ray counterparts (GWR + X) and both (GWR + Opt + X). Similarly the top right panel shows the population that is in the direct 
impact phase of accretion in grey and we subdivide that population in GWR and X-ray sources. The bottom two panels show (again in grey) the populations 
that are detectable in the optical band (left panel) and the X-ray band (right panel). The distribution of sources detectable in both the optical and X-ray band is 
shown by the hatched histogram in both lower panels (Opt + X). See also Tablel3l 
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Table 3. The number of short period (P < 1500 s) AM CVn systems with 
different properties. The number of systems with each property that is in the 
direct impact phase is given in the second column. See also Fig.|9|for the 
distribution of these systems over orbital period. 



rate is very sensitive to small changes in assumptions in binary 
population synthesis calculations. One has to keep in mind that 
a lot of these assumptions are badly constrained, which means 
that the number of p rogenitor systems is rather uncertain (e.g. 
iNelemans et al 1200 id) . On top of that, there are uncertainties con- 
nected with the formation of AM CVn syste ms from their imme- 
diate progenitors as discussed in Sect. l2.ll and NPVY, making the 
AM CVn birth rates uncertain by at least an order of magnitude. 
The models we use for the X-ray and optical emission are very 



simplified as well. The interstellar absorption model might very 
well underestimate the X-ray absorption for the systems in the 
Galactic centre, which comprise a large fraction of the systems that 
are detectable only in X-rays (pluses in Fig.0. We furthermore re- 
call that a number of (possibly important) effects like irradiation 
and compressional heating are not taken into account in our simple 
emission models. 

For comparison we also computed the results for X-ray de- 
tectors sensitive in the range 0.1 - 15 keV (similar to the Chan- 
dra and XMM bands). Because most of our spectra are rather soft, 
the flux (T) in this band is generally not much larger than in the 
ROSAT band: 80 per cent has Jxmm/^rosat < 1-5; 96 per cent 
has J-xmm/^rosat < 3. However, the improved sensitivity of 
Ch andra and X MM obviously will have an impact. For instance 
the Muno et al. (2003) Chandra observations of the Galactic Cen- 
tre have a completeness limit of 3 10 -15 erg cm -2 s _1 , almost 
two orders of magnitude deeper than our assumed ROSAT limit. 
The number of X-ray sources in the 0.1-15 keV band detectable 
down to 10~ 14 erg cm" 2 s _1 is 644 and 1085 down to 10~ 15 
erg cm~ 2 s" 1 . In the Chandra mosaic image of the Galactic centre 
( Wang 'et all2002l) . roughly down to 10~ 14 erg cm -2 s _1 there are 
~1000 point sources, while our model gives 16 X-ray systems in 
this region. In the tiny Muno et al. 1 2003) region, where they iden- 
tify more than 2000 sources, we predict 3 AM CVn systems with 
X-ray flux above 10~ 15 erg cm~ 2 s -1 , 
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As briefly discussed in Sect. [2] we did not include AM CVn 
systems formed like CVs in our models. The reason is that for 
the production of the short period systems, which we consider in 
this paper, this channel is not effective. An evolutionary track for a 
donor with an initial mass of 0.96 M s and a hyd rogen-exhausted 
core ~ 0.05 Mq (sequence nr. 4 in iFedorova & Ergmalll989Tl 
reaches a minimum period of 8.4 minutes and spends 3.8 10 7 yr 
in the 25 - 8.4 minute period range (A.V. Fedorova, private com- 
munication). Taking this time as typical, and assuming a birth rate 
of 6.4 10~ 5 yr -1 for systems evolving to P < 25 minutes, as in 
the standard model of Podsiadlo wski et ail (2003), we expect 2400 
post-CV ultra-compact systems in the Galaxy, compared to the 
140,000 from the channels discussed here. However, if these chan- 
nels are ineffective, e.g., due to the instability of the mass transfer 
and/or destruction of accretors by edge-lit detonations l NPVY|), the 
AM CVn systems from CVs may become more important even at 
short periods. 



tion for the donor stars in these systems. LISA is expected to detect 
a similar number of detached double white dwarfs, for many of 
which the frequency change can be determined, yielding a direct 
measurement of the distances to these sources. Maybe for a few the 
second derivative of the frequency can be measured, which would 
give a direct test whether the frequency change is purely due to 
GWR, or whether tidal effects play a role. 

Our results indicate that many short-period AM CVn systems 
may be detectable in more than one way. A total of 106 systems is 
detectable in the optical and X-ray bands. For the systems that are 
detectable with LISA, about 300 (3 per cent) have optical and/or X- 
ray counterparts, giving the opportunity to independently determine 
the position and orbital period. This leads to a better determination 
of the other parameters that can be derived from the GWR signal 
(i.e. inclination and masses) and would lead to unprecedented ac- 
curacy in the study of these extreme binaries. 



7 CONCLUSION 

We discussed a model for the Galactic population of short-period 
AM CVn systems and the uncertainties in this model. We find a 
total population of AM CVn systems in the Galaxy with periods 
shorter than 1500 s of ~ 140,000 systems. We estimate that the fi- 
nal uncertainty on this number is at least a factor of 10. The model 
presented here gives an optimistic estimate of the total number of 
systems, i.e. all quoted numbers could be (substantially) smaller. 
We calculated the optical, X-ray and gravitational wave signals that 
could be detectable, using simple emission models and selection ef- 
fects, and discussed the overlap between the subpopulations which 
may be discovered by different detection methods. 

In this optimistic model we find that there are a few systems 
expected that are in the direct impact phase and are detectable as 
optical and X-ray sources, where the optical emission comes from 
the cooling of the (possibly heated) white dwarfs. Thus if RX J0806 
and V407 Vul turn out to be short period AM CVn systems, that 
would argue for an optimistic model as in a pessimistic model 
there woul d be no such system s that are observable. Stro hmaveil 
(2002|) and lHakala et"ail 12003) have reported a shortening of the 
periods in V407 Vul and RX J0826, arguing against this interpre- 
tation despite the fact that these stars might be helium-donor or 
CV sys tems before the period minimum (see also Wpudl& Warner 
(2003a) who argue that there is not sufficient evidence for a short- 
ening of the period in RX J0806). 

At the observed periods of RX J0806, V407 Vul and ES Cet, 
the fraction of systems that is in the direct impact phase is 80, 40 
and 30 per cent respectively. The total number of direct impact sys- 
tems in the Galaxy is 2680. We further find that there are 220 (very) 
short period AM CVn systems that are only detectable in the X-ray 
band, some with orbital periods as short as 3 minutes. 

We conclude that the future gravitational wave detector LISA 
is expected to discover thousands of short-period AM CVn sys- 
tems, about ten per cent of the total number of short-period systems 
in the Galaxy. The number of systems with periods below 1500 s 
that is detectable in the optical band is 1336 (about one per cent of 
the total). Only 326 systems (0.25 per cent) are detectable in the 
ROSAT X-ray band. For the typical sensitivities of Chandra and 
XMM, the number of detectable sources over the whole sky is > 
1000, still only a tiny fraction of all X-ray sources. 

For many of the LISA sources the change in frequency can be 
measured, which would provide a test of the mass - radius rela- 
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Figure Al. Cooling curves of white dwarfs used in this paper. The curves 
are a rough translation of Fig. 10 o flHansenlfl999l) . 
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APPENDIX A: WHITE DWARF COOLING 



We mo del the cooling of white dwarf based on the curves o fH ansen 
( 1999). We abandoned the previous use of the D riebe et alJ 1 1998) 
models, since their models with thick hydrogen envelopes cool 
slower than observations of detached double white dwarfs suggest 
and their thick hydrogen envelopes actually may be lost in hydro- 
gen shell flashes (e.g. lNelemans et all2001ct) . 

Our simple models have carbon-oxygen white dwarfs that cool 
faster than helium white dwarfs, but within the groups the most 
massive objects are brighter at given age. We model the curves as 
(see Fig. IAU : 

logL = L max - 1.33 log(£/10 6 yr), (Al) 
where L max is a linear fit given by 

Lmax = l-(0.9- M WD ) for A/wd > 0.5M Q , (A2) 

and 

L max = 1.4 - 1.33(0.45 - Mwd) for M W d < 0.5M Q (A3) 

(mass and luminosity in solar units). For white dwarf masses below 
0.5 Mq the luminosity is constrained to be below log L/Lq = 
—0.5, for more massive white dwarfs below log L/Lq — 2. 
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